Fluctuational phase-flip transitions in parametrically pumped oscillators 
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We analyze the rates of noise-induced transitions between period-two attractors. The model investigated 
is an underdamped oscillator parametrically driven by a field at nearly twice the oscillator eigenfrequency. 
The activation energy of the transitions is analyzed as a function of the frequency detuning and field 
amplitude scaled by the damping and nonlinearity parameters of the oscillator. The parameter ranges 
where the system is bi- and tristable are investigated. Explicit results are obtained in the limit of small 
damping (or strong driving), and near bifurcation points. 
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I. INTRODUCTION 

Nonlinear systems driven by a sufficiently strong pe- 
riodic field often display period doubling The two 
emerging stable periodic states are the same: the only 
difference between them is that they are shifted in time 
by the period of the force ^-k jtop. This feature is a con- 
sequence of the symmetry with respect to translation in 
time by 2'it/ljJp, and it has attracted much attention to 
such systems as elements of digital computers. In the 
presence of noise, there occur fluctuational transitions 
between the period- two attractors, which correspond to 
phase slip of the system by tt. Such transitions have 
been observed recently for electrons oscillating in a Pen- 
ning trap 1^. They have been also investigated in an 
analog electronic circuit in the context of stochastic 
resonance 

Motivated by these observations, in the present paper 
we develop a theory of escape rates from period-two at- 
tractors. The analysis is done for the simplest generic 
model that displays period doubling: an underdamped 
oscillator parametrically driven by a force at nearly twice 
the oscillator eigenfrequency ujo ||]. This model applies, 
in particular, to axial vibrations of an electron in a Pen- 
ning trap . Much work on a parametrically excited os- 
cillator has been done in the context of squeezed states of 
light, cf. H]. We analyze escape due to classical fluctua- 
tions, which were substantial for the systems investigated 
in|i. 

Parametrically excited oscillator is an example of a sys- 
tem away from thermal equilibrium. Such systems usu- 
ally lack detailed balance 0, they are not characterized 
by free energy, and escape rates depend on the system 
dynamics and the noise that gives rise to fluctuations 
in the system. In the important and quite general case 
where the noise is Gaussian, there has been developed a 
technique which reduces the problem of calculating es- 
cape rates to a variational problem [^. The solution of 
this problem describes the optimal path along which the 
fluctuating system is most likely to move when it escapes. 
Such path is often called the most probable escape path 



(MPEP) . The value of the variational functional gives 
the exponent in the expression for the escape rate. 

An advantageous feature of the problem of period dou- 
bling in an underdamped oscillator is that the period 
doubling occurs for comparatively small amplitudes F of 
the driving force, where the nonlinearity of the oscillator 
is still small: the anharmonic part of the potential energy 
is much less than the quadratic term Lo^q^ /2. In this case, 
the quantities of interest are the amplitude and phase of 
the vibrations at the frequency uf/2 « loq. They vary 
only a little over the time ~ w^^. The corresponding 
dynamics is affected by Fourier components of the noise 
within a narrow band centered at LOpj^. Essentially, this 
means that, in the analysis of the dynamics of slow vari- 
ables, the noise may be assumed white. A similar situa- 
tion arises pO| in the problem of transitions between the 
stable states of forced vibrations of a resonantly driven 
underdamped oscillator. 

Below in Sec. II we provide the results on the the phase 
portrait of the oscillator in the rotating frame, derive the 
properties of noise for slow variables, and, for low noise 
intensities, formulate the variational problem for the ac- 
tivation energy of escape. In Sec. Ill we provide explicit 
expressions for the escape rates in the following parame- 
ter ranges: vicinities of the bifurcation points where pe- 
riod two attractors emerge and where there emerge two 
unstable period two states, and the range of compara- 
tively strong driving where the motion in slow variables 
is underdamped. In Sec. IV numerical results for the gen- 
eral case are discussed and compared with the results for 
the limiting cases. Sec. V contains concluding remarks. 



II. ASYMPTOTIC EXPRESSION FOR THE 
ESCAPE RATES 

A. Phase portrait in slow variables 

To set the scene, we will first discuss the phase por- 
trait of a parametrically driven underdamped oscillator. 



A simple phcnomenological equation of motion is of the 
form 



dt^ ^ dt 
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qF cos oj pt ^ £,{t) (1) 



Here, F is the friction coefficient, 7 is the nonhnearity 
parameter, F is the amphtude of the regular force, and 
^{t) is a zero-mean noise, {^{t)) = 0. The effect of cubic 
nonlinearity in the Hamiltonian of the oscillator can be 
taken into account by renormalizing the parameter 7 (cf. 
1^ ) . The results can be also easily generalized to the case 
where, for special reasons |^], 7 is numerically small, and 
it is necessary to allow for the term cx q^ in the equation 
of motion. 

We consider resonant driving, so that the oscillator 
eigenfrequency loq is close to lofI2, 



r, |2a;o - wfI < <^vi- 



(2) 



In this case it is convenient (cf. g) to analyze the os- 
cillator motion in the rotating frame. Respectively, we 
change to slow dimensionless time t = and slow di- 
mensionless variables qi,q2- 



dq 
di 
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(3) 



Following the standard procedure of the method of 
averaging (cf. and neglecting fast oscillating terms 
which depend on the oscillator amplitude and contain a 
factor exp(±i7iCt;i?i/2) with n ^ 0, one obtains the equa- 
tions of motion for q\ , (72 in the form 



qi ^ -J- = -qi + ^ ^ Ci(t/F), t = Tt, 

dr oq2 

■ - ^^92 dg 

92 = = -92 - h 6(t/F) 

dr oqi 
where ^1,2 (t/F) are the random forces, and 

ql) 
-ql) 



(4) 



5(91,92)= 2(^1 



{ql 



92jsgn7 



(5) 



7tC(9| 



In what follows we assume that 7 > 0; the case 7 < can 
be described by replacing Vl —17, qi —qi {i — 1, 2), 
and by changing the phase of the field in Eq. (^ by tt. 

Except for the random force, the motion of the oscil- 
lator as described by Eqs. (Q) is characterized by two 
dimensionless parameters; the scaled frequency detuning 
fl and the scaled field C, 



n=[{ujF/2)-uja]/T, C,=F/2ujpT. 



For C < 1 or for < -(C^ - 1)^/^ the osciUator ^ in the 
absence of noise has only one stable state, qi = q2 = 0: 
parametric oscillations are not excited. The value C = 1 
gives the threshold field amplitude Fth = 2u!p which is 
necessary for exciting period-two vibrations. The phase 
portrait of the oscillator in variables qi , 92 in the range 
where the oscillations are excited is shown in Fig. 1. 

For > 1, with the increasing Q there first occurs a 
pitchfork bifurcation for fl = — (C^ — 1)^/^. Here, the sta- 
ble state 9i = 92 = becomes unstable and there emerge 
( 1 2) • 

two stable states q^^' which are symmetric with respect 
to qi — q2 = 0, see Fig. la. These states correspond to 
stable period two vibrations shifted in phase by tt. The 
vibrational amplitude a = \qst \ increases monotonically 
with the increasing SI (see inset in Fig. 1). The results of 
the asymptotic analysis based on Eqs. (^) apply for not 
too large amplitudes a, 

Ta^ ,\2ujp — uiqIo^ ojp. (7) 

As fl goes through the value (C^ — 1)^/^ there occurs 

the second pitchfork bifurcation of the state 91 = 92 = 0: 

this state becomes stable again, and there emerge two 
(1 21 

unstable states , which are also symmetric with 

respect to the origin and correspond to unstable pe- 
riod two vibrations of the oscillator. As it is seen from 
Fig. lb, the phase plane of the oscillator is such that, for 
SI > (C^ — ly^^ the domains of attraction to the stable 
states q^P , c{^^^ are separated by the domain of attraction 
to the stable state q = 0. 



B. Noise in the equations for slow variables 

In many cases of physical interest, the noise ^(i) driv- 
ing the oscillator is Gaussian. This noise may originate 
from the coupling to a thermal bath (which also gives 
rise to the friction force in Eq. (|l|)), or it may be due to 
an external nonthermal source. A zero-mean Gaussian 
noise is characterized by its power spectrum 

/•OO 



(8) 

For a stationary noise the power spectrum (^) is indepen- 
dent of time. In what follows we assume that the function 
(j){uj) is smooth near the oscillator eigenfrequency loq- 

The noises £,i,2(t) in Eq. (^ are nonstationary, gener- 
ally speaking. One obtains from Eqs. (^, (||), (^ the 
following expressions for their power spectra: 



3|7l 
8cj3F 
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(6) ^u[S.i{t + t'),^2{t)V 



X [1 - (-1)* cosujpt] {i ^ 1, 2), (9) 

3|7l J 1 \ . 
— j-^(p [ (u — T^^J-T ) smujpt. 



It follows from Eq. (|^) that the spectra of the diagonal 
correlators {£,i{t)(,i{t')) have both time-independent and 
fast-oscillating in time components, whereas the power 
spectrum of the cross-correlator of , ^2 is fast oscillat- 
ing in time. Therefore, in the analysis of the effect of the 
noise on the slowly varying functions qi.2, in the spirit 
of the averaging method, one can assume that the noise 
components ^lit), ^2{t) are asymptotically independent 
of each other; one can also leave out the terms oc cos LOpt 
in the power spectra of the diagonal correlators (cf. ||ic| ] 
where the analysis was done for a microscopic model of 
noise resulting from coupling to a bath; rigorous math- 
ematical results on the method of averaging in noise- 
driven systems are discussed in [|lT|). 

The dynamics of the slow variables qi^2 is character- 
ized by the time scales ~ l/F, 1/\2u!p — ujq\. If the power 
spectrum 4>{uj) varies only slightly in the whole range 
\uj — ujp/2\ < r, 1 2ti;i? — I (and this range does not corre- 
spond to a deep minimum of then one can assume 
that the random functions ^1 (t) , ^2 (t) are asymptotically 
independent zero-mean Gaussian white noises, 



of a variational problem. This problem is of the form (cf. 
00): 



Ur/m{r'/r))^D~S{T~ 
D = (3|7|/8a;|,r2)^(wj./2) 



1,2 



(10) 



(the function 5{x) is of the type of a (5-function: it is large 
in the domain ^ 1, and its integral is equal to 1). 

The quantity D gives the characteristic intensity of the 
noise in the equations of motion for slow variables qi , q2 
(^. In particular, if the oscillator is coupled to a thermal 
bath with the correlation time much smaller than I/wq, 
so that it performs a "truly" Brownian motion, and the 
random force £^{t) in Eq. (|l|) is ^-correlated, with intensity 
4rfcT, we have D = 3|7|fcr/2w|,r. We note, however, 
that the dynamics of slow variables can be described as 
Brownian motion even where this description does not 
apply to the motion of the initial oscillator, i.e. where the 
correlation time of the thermal bath is ^ l/wo- In this 
latter case, the friction force in Eq. (^ is also retarded 
(cf. @). 



C. Activation energy of escape 

If the dimensionless noise intensity D is small, for most 
of the time the oscillator fluctuates in a small vicinity 

in) 

of one or the other stable state q^j . Only occasionally 
there occurs a large fluctuation which results in a tran- 
sition to another stable state. The probability Wn of 
such fluctuation is exponentially small, and its depen- 
dence on the noise intensity is given by the activation 
law, Wn oc exp(— ^n/Z?) (see |^ for a review). In fact, 
to logarithmic accuracy Wn is determined by the prob- 
ability density of the least improbable realization of the 
force $,{t /T) which results in the corresponding transi- 
tion. Therefore the quantity S'„ is given by the solution 



Wn = CeM-Sn/D), Sn = min5[q(T)], (11) 

/oo 1 
drL(q,q), i(q, q) - - [q - K(q)]^ 
-00 ^ 

q(t) — > q^"^ for t ^ — cx), q(t) q„ for t — > 00, 

where the components of the vector K are given by the 
rhs parts of the equations of motion (H), 



^i(q) 



dg_ 
dq2 ' 



K2{^) 
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dg 



(12) 



(in what follows we set n = 1,2 for the stable states of 
period two vibrations, and n = for the stationary state 
q = where it is stable) . 

The solution of the variational problem ([l^) q(i) 
describes the optimal, or most probable escape path 
(MPEP) from the nth stable state. This path is 
instanton-like (see 0|): it starts at the stable state q^""* 
for t — > —00, and for t ^ 00 it approaches the un- 
stable state q„ on the boundary of the domain of at- 
traction to q^"^ (having reached the boundary, the sys- 
tem makes a transition to another stable state with a 
probability ~ 1/2). It follows from Fig. 1 that, for 

_(^2 _ ]^-)l/2 < < (C^ - 1)1/2, ^jjg Qj^ly 

ble states are period two attractors, escape from one of 
them means a transition to the other. For > (C^ — 1)^/^ 
escape from one of period two attractors means a tran- 
sition to the state where period two vibrations are not 
excited. From this state, the system makes fluctuational 
transitions into one or the other state of period two vi- 
brations. 

The most probable realization of the force is related 
to the optimal fluctuational path via Eq. (||), ^(r/F) = 
q — K. Optimal fluctuational paths are physically real, 
they have been observed in the experiment fl^ . 



III. ACTIVATION ENERGY OF ESCAPE IN THE 
VICINITIES OF THE BIFURCATION POINTS 



The activation energies Sm as defined by (|ll|), depend 
on two dimensionless parameters of the driven oscillator: 
the scaled field strength C, and the scaled frequency de- 
tuning (^). In the general case, Sn may be calculated 
numerically (see Sec. VI). Explicit expressions for Sn can 
be obtained in several limiting cases. 

In the range of comparatively small nonlinearity (|^), 
the oscillator may experience two bifurcations with the 
varying field frequency, as seen from Fig. 1. In the vicin- 
ity of a bifurcation point, one of the motions of the system 
near the emerging stable state(s) becomes slow: there 
arises a "soft mode" [Q. Correspondingly, fiuctuations 
near a bifurcation point have universal features |15|. 



Near a bifurcation point one can either solve the vari- 
ational problem (^T]) explicitly or reduce the system of 
equations of motion to one first order equation. From 
this equation one can find not only the exponent, but also 
the prefactor in the expression for the escape rate. 

The equation for a slow variable Q can be derived from 
@ by appropriately rotating the coordinates: 

Q = qi cos (3 + q2 sin f3, P = —qism(3 + (72 cos /3 (13) 
tan2/3 = ng = t(C^ - 1)^^^ 

(fls are the bifurcation values of the dimensionless fre- 
quency detuning, see Fig. Ic). For jfi — fisj the 
dimensionless relaxation time of the variable Q goes to 
infinity, whereas that of P is 1/2, and therefore P{t) 
follows the slow variable Q{t) adiabatically [|l5|. The 
fluctuations in P can be neglected compared to the fluc- 
tuations in Q. Using the adiabatic solution for P (i.e., 
neglecting P and the noise term in the equation for P) 
we obtain the following equation for Q: 

Q = -dU/dQ + S(r), (S(t)S(t')) - D~S{t - r') (14) 
■l_ „ 1 



[/(g) = 



-(r!-r!B)g2--cV 



, D.b\^-^b\ < 1. 



It is seen from Eq. (|T^) that, near the bifurcation point 
^Ib = — (C^ ~ l)^^^i the system has one stable state for 
il < and two symmetrical stable states for $7 > il^ 
(cf. Fig. 1). The escape rates from these symmetrical 
states are the same and are given by the Kramers ex- 
pression p3] 



= (V2/7r)|r!s|(r! - ^b) exp(-5„/i?), (15) 

Sn = 2 [{/(Q(")) - c/(q("))] = \nBm - ^Bfi^e 

(Q(") and Q(") are the values of the coordinate Q in the 
stable and unstable states). 

The rate (|l^) is the rate at which there occur phase 
slip transitions between the period two stable states for 
n close to ^B = -iC - 1)^/^- Eq. (|l|) applies for 
exp(— S'n/D) <C 1. The activation energy Sn is quadratic 
in the distance |f2 — ^b\ to the bifurcation point. For a 
given |rj — fJsl the dependence of Sn on the dimensionless 
field ( is given by the factor {(^ — 1)^/ j^^. 

For < 17 - (C^ - 1)^/^ < 1, Eq. (|l|) describes the 
rate of transitions from the stable state q^i = to each 
of the period two stable states. The transitions occur via 
the appropriate unstable period two state (the one on 
the boundary between the domain of attractions to the 
stable period two state and the stable state q = 0). 



IV. ACTIVATION ENERGIES OF ESCAPE IN 
THE SMALL-DAMPING LIMIT 



A. Motion in the absence of dissipation 

Of special interest is the case where the scaled field 
amplitude C is large enough so that the dissipation terms 
— gi,— 52 in the rhs of the equations for slow variables 
(^) are comparatively small. In the neglect of these terms 
and the random force, Eqs. (^) describe conservative mo- 
tion of a particle with the coordinate qi and momentum 
q2, and with the Hamiltonian function g{qi, 92) @- This 
particle moves along closed trajectories shown in Fig. 2 
for the auxiliary variables X,Y: 



X = q,/e/\ Y^q^/C 

dX dG dY dG 



1/2 



dX' 



G{X,Y) 



-- (r. (16) 

[c'/'xx'^'y)^Io,-i)x' 



df dY' df 



C'9 



-lil. + l)Y'-\iX' + Yr,f^-^. 



The conservative motion (|T^) depends only on one pa- 
rameter, fi = ri/C, which characterizes the interrelation 
between the frequency detuning and the field strength. 
The shape of the effective energy G{X, Y) for two differ- 
ent values of fi is shown in Fig. 3. The trajectories in 
Fig. 2 are just the cross-sections of the surface G{X, Y) 
by planes G =const. The extrema of the surface G(X, Y) 
are the fixed points of the system. 

For /i < — 1 the surface G{X,Y) has one extremum 
placed a,t X — Y = 0, which corresponds, with dissi- 
pation taken into account, to the stable state with no 
period two vibrations excited. For — 1 < /i < 1 the func- 
tion G{X, Y) has two maxima at X = 0, Y ^ ±{ij.+ iy/^ 
(they correspond to two period two attractors) and a sad- 
dle point a,t X = Y — 0. For /i > 1, in addition to the 
above maxima the function G{X, Y) has a minimum at 
X = Y = (the maxima and the minimum correspond to 
the stable states of the oscillator) , and two saddle points 
at X — ±{fj, — 1)^/^, y = 0. The extreme values of G are 
given by the expressions 



Gl.2 = ^(/i+l)' 



Gu = for /i < 1; 



(17) 



Go — 0, Gu 



-(^-1)2 for^ > 1. 



We note that for /i > 1 the trajectories surrounding the 
states at X = 0, Y — ±(/i-|-l)^/^ in Fig. 2 become horse- 
shoe like for large enough Gi — G, and also that for < 
G < Gu there are coexisting "internal" and "external" 
trajectories. 



B. Escape rates 

The effect of small dissipation in Eqs. (^) is to trans- 
form the closed trajectories in Fig. 2 into small-step spi- 
rals which wind down to the corresponding stable states 



(cf. Fig. 1). The motion can be described in a standard 
way in terms of slow drift over the energy g towards the 
stable state (cf. |^]). 

The random force which drives the system away from 
the stable state should "beat" this drift. It would be ex- 
pected that the optimal fluctuational path corresponds 
to energy diffusion away from the stable state. A solu- 
tion of the variational problem (^) for small dissipation 
was obtained in [p^ for a different form of the function 
5(91? 92)- Alternatively, one can use the approach first 
suggested by Kramers in the analysis of escape of un- 
derdamped thermal equilibrium systems. This approach 
is equivalent to deriving from Eqs. an equation for g 
and averaging the dissipative term and the noise inten- 
sity in this equation over the period of vibrations with 
a given g in the absence of dissipation and noise. The 
dissipative term (the dissipation rate of g) is then given 
by the expression —[qi{dg/dqi) + q2{dg/dq2)], with an 
accuracy to the correction oc D (here, overline means 
averaging over the vibration period). The diffusion co- 
efficient for g is given by D[dg / dqiY + {dg/dq2)^]- The 
resulting first order equation for g can be solved to give 
the following expression for the activation energy: 



Sn = 2C r^G^^, M{G) 

JCr, 



N{G) = - 



A(G) 



N{G) J jA^a) 

dXdYV^G{X, Y). 



dXdV, 



Using the explicit form of G(X, Y) one obtains 



NiG) 



[ [ dXdY [fi - 2{X^ + Y^)] 

J J A{G) 



(18) 



(19) 



The values of G„ {n =1,2) and G„ are the extreme val- 
ues of G(X, Y) as given by Eq. (|l^). The double inte- 
grals in (|l^) are taken over the areas A{G) limited by 
the trajectories G(X, y) = G in Fig. 2 which surround 
the nth center (the expressions for A/(G), N{G) were ob- 
tained from the expressions for the drift and diffusion 
coefficients for g using the Stocks theorem, with account 
taken of Eqs.(|l|), as it was done in [ M). 

The expressions for M{G),N{G) (^ can be fur- 
ther simplified by changing to polar coordinates X = 
R cos ip, Y — Rsimp. Solving Eq. (16) for in terms of 
G, (p and integrating over i?^ one then obtains that, in 
the problem of the escape from the period two attractors 
(n = 1, 2 in (|l|), i.e., G > G„, cf. Fig. 3), 

Af(G) = I d^/(G,(p), (20) 
N{G) = - j dip i^i -2cos2ip) f{G,cp), 

/(G,(p) = [(m- cos 2(^)2 -4G]'^' (G>G„). 

The limits of the integrals over ip are determined from 
the condition that /(G, (p) be real and /i — cos2<y9 > 0. 



In the problem of escape from the stable state at q = 
for ^ > 1 (G < G„, see Fig. 3) 

Af(G) = J dpR^{G,p), (21) 

N{G) = J d>pRHG,^)[fi-R\G,p)] (G<G„) 

R^iG,p) = n-cos2p- [(^ - cos 2^)2 -4G]^^^ . 
(here, the integral is taken from to tt). 

C. Explicit limiting expressions for the escape rates 

The expressions for the activation energies are simpli- 
fied near bifurcation points (yet not too close to the bifur- 
cation points, so that the damping of the vibrations with 
a given g ( p^ be small compared to their frequency). 
These bifurcation points are /i = ±1. 

As the increasing /i goes through the value /i = — 1, 
the maximum of the function G(X, Y) &% X = Y = Q 
becomes a saddle point from which there are split off two 
maxima of G corresponding to period two attractors (see 
Fig. 3). With the further increase in /i, for /i = 1 the 
point X = Y = Q becomes a minimum of G{X, Y) from 
which there are split off two saddles of G. Escape from 
the emerging stable states is determined by small-radius 
orbits G(X,Y) = G. For such orbits M{G)/N{G) w l//i 
in (fsh, and therefore 



Si ^82^ dp + 1)V2, < + K 1, 
^o«C(m-1)V2, 0</i-l«l. 



(22) 



In the limit of large fj,, the activation energy for escape 
from period two attractors (|l^) is determined by orbits 
with G„ = (/i - 1)2/4 < G < Gi,2 = (m + 1)V4. These 
orbits have a shape of narrow arcs on the {X, y)-plane. 
It is seen from ^ that, for such orbits M{G)/N{G) « 
l//i, and 



81^82^ 2C, /i > 1- 
One can also show from (|2^) that 

80 ^ CfJ-, M > 1- 



(23) 



(24) 



In the general case of arbitrary values of the parame- 
ter fi — n/( the activation energies 81 — 82 , 80 could be 
found by evaluating the integrals (p^), (pO|), ( pl| ) numer- 
ically. The results are shown in Fig. 4. 

It is seen from Fig. 4 that 5*1 = 82 is quadratic in 
{lJ,+ l) for small 11+ I, and monotonically increases with 
the increasing fi. For large fi, 81,2 saturates at « 2C . 
The activation energy 80 is quadratic in /i — 1 for small 
/i — 1, and then monotonically increases with the increas- 
ing fi. We note that 81 > 80 for fi < 6.5. Respectively, 
for such /i the escape rates from the period two attractors 



are exponentially smaller than from the state where the 
vibrations are not excited. As a consequence, the station- 
ary population of the period two attractors wi — W2 is 
exponentiaUy larger than the population of the steady 
state: 

wi^W2 = {Wq/2Wi)wq, wi/wo oc exp[(S'i - So)/D]. 

(25) 

For larger frequency detuning {Q = /i^), the steady 
state becomes more populated, in agreement with the 
intuitive physical argument that, as the field is detuned 
further away from the resonance, it is less likely that the 
period two vibrations will be excited, for the same field 
intensity. 



V. EFFECTS OF THE SIXTH-ORDER 
ANHARMONIC TERM 

In the experiments ^ , because of the structure of the 
electrostatic field in the trap for an oscillating electron, 
the 6th order anharmonic term in the Hamiltonian of the 
electron vibrations could be relatively large (in fact, the 
4th order term could be relatively small), while higher- 
order terms remain much smaller than both the 4th and 
6th order terms [ p^ . The advantage of suppressing the 
4th order term is that the amplitude of the period-two 
vibrations becomes larger. When the 6th order anhar- 
monicity is taken into account, the equation of motion 
takes the form 



2r 



dq 
'di 



-ujlq- 



^q^ + Xq^ + qF cos ujFt = ^{t) 

(26) 



For the model (|26D, the equations of motion in the 
rotating frame are again of the form (|^), but now the 
function 3(91,92) is given by the expression 



5(91, 92) = ^^{ql + gf) - ^{ql + g|)^sgn7 



(27) 



P / 2 , 
5XF 



p 
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We have neglected the renormalization (c>c 7^) of the non- 
linearity parameter parameter A (this renormalization is 
substantial when 7 is not small, in which case the role 
of the 6th order anharmonicity is insignificant, or the 
whole approximation of small-amplitude vibrations does 
not apply). The dimensionless parameter p characterizes 
the "strength" of the 6th order nonlinearity. 

For fi very close to the bifurcation values t(C^ ~ l)^^^j 
escape from the small-amplitude stable state(s) is deter- 
mined by the motion with small q\ + gf. Clearly, this 



motion is determined by the lowest- order anharmonic- 
ity, and therefore neither this motion nor the bifurcation 
values of the parameters are affected by the higher-order 
nonlinear terms. 

In what follows we will investigate the effect of the 



term oc A in (26) on the escape rates from the period 



two attractors, and also from the zero amplitude stable 
state for il > (^^ — 1)^^^7 far from the bifurcation points. 
We will analyze the case of weak damping, which is of 
utmost interest for the experiment [0. 

In the neglect of dissipation and fluctuations, the mo- 
tion of the oscillator ( |26| ) in the rotating frame can be 
described by Eqs. (|l6|), but when the 6th order nonlin- 
earity is taken into account, the effective energy G{X, Y) 
is of the form 

\{p-l)X^ + 



G{X,Y) 



2' 

y2Vl 



p{X- 



1 

2 

+ r2)3 



(A^ + l)^' 



p 



(28) 



For p > (i.e., for 7A > in Eq. (26|)), the phase 
portrait of the conservative motion (|l6), ( p8| ) remains 
qualitatively the same as that shown in Figs. 2a,b for 
p = 0, as does also the topological structure of the sur- 
face G{X^ Y) in Fig. 3. The centers which correspond 
to the period two attractors lie on the F-axis, Xx.2 = 0, 
they correspond to the maxima of the function G{X, Y), 
whereas the hyperbolic points which correspond to the 
unstable period two states lie on the X-axis, l^ni,2 — 0; 
they correspond to the saddles of G{X,Y). 

The general expression for the activation energy of es- 
cape Sn is given by Eq. ([l8|), but now the functions M{G) 
and N{G) have to be calculated with account taken of 
the 6th order terms in G{X^ Y) (^). In particular, 

N{G)^ 1 1 dXdY [/i-2(A:2 + r2) 
J Ja{g) 

- 3p{X^ + y2)2] (29) 

The functions M{G),N{G) can be written as single 
integrals over the polar angle (cf. Sec. IV B). In the 
problem of escape from period two attractors, we obtain 
that, similar to (BC 



MiG) = \j [R\{G,v) - RliG,^)] 



(30) 



N{G)^- J dip [Z+(G,^)-Z_(G,^)], 

Z±{G,^) ^ pRliG,^) ^ RiiG,^) ~ pR%{G,^) 

Here, R±{G,Lp) are the radii of the trajectories 
G(A", Y) = G which surround the centers correspond- 
ing to the period two attractors. They are given by the 
real roots of the equation 

lpRl + ]Ri~l-ip~cos2ip)Rl + G = 0, (31) 
6 4 2 

Gi,2 > G > for /i< 1, Gi^2 > G > G„ for /i > 1. 



(Gi — G2 are the values of G in the period two attractors, 
see Eq. (34)). The limits of the integrals over ip in Eq. ( |30| ) 
are determined from the condition that Eq. ( pi] ) has two 
real roots, i?_ < R+. 

For the case of escape from the attractor with zero vi- 
bration amplitude (q=0) the hmctions M{G) and N{G) 
have the form 



1 /■^'^ 



(32) 



N{G)^- dipZ_iG,ip), M>1, >G>0 
^ Jo 



where the function Z-{G, ip) is defined in Eq. ( |3C| ) (here, 
G„ is the value of G in the saddle points, see Eq. (34); 
in the range < G < Gu Eq. ( |3l| ) has only one real root 
which determines the function Z-). 



1. Strong sixth order nonlinearity 

The activation energies (|3^ ) as functions of the dis- 
tance /i^ — 1 to the bifurcation points display an interest- 
ing behavior for large sixth-order nonlinearity, p 3> 1 , in 
the range where /.t^ — 1 <C 1 but p{^'^ — 1) ^ 1. It follows 
from (|3^), (HJ) that in this range 



^i,2~:7C(m + 1)^^V"'/', M+Kl, p(/i + !)>!, (35) 



^0 ~ ^Cif^ - l)3/'p-l/^ M - 1 « 1, p(m - 1) » 1- 

The dependence on the distance to the bifurcation point 
/i^ — 1 as given by Eq. (35) is described by the power law 
with the exponent 3/2. In contrast, for weak sixth order 
nonlinearity ( p2|) the exponent is equal to 2. Eqs. (|3^), 
( p^ ) describe both limiting behaviors and the crossover 
from one of them to the other. 



A. Activation energies of escape near bifurcation 
points 



B. Strong sixth order nonlinearity: general case 



For underdamped systems where the sixth order an- 
harmonicity, the analysis of the activation energies of es- 
cape for parameter values close to the bifurcation points 
is similar to that in Sec. IV C. In the range < 1 <C 1 
the centers (-'^^1,2, ^1,2)) which correspond to the period 
two attractors are close to the saddle point at the origin. 
In this case R± <C 1 in Eqs. (|3l]), and the terms with 
_|_ y2 jj-^ ( p9| ) can be neglected. These terms can 
be also neglected in the problem of escape from the sta- 
ble state X = y = for < /i -f 1 ^ 1, as in this case 
the hyperbolic points (X„i_2, ^^111,2) are close to the sta- 
ble state. Therefore, in both cases, M{G)/N(G) « l/^i 
(cf. Sec. IV C), and it follows from Eq. (|l|) that the 
activation energies of escape 



S-i = 52 « 2CGi, 0<M+1<1, 
Sq w 2CG„ < - 1< 1 



(33) 



Here, we have taken into account that G(0, 0) = 0. The 
explicit expressions for the effective energies Gi = G2 
(the maxima of G{X,Y)) and G„ = Gui — Gu2 (the 
saddles of G{X,Y)) are of the form 

r Kl + 6p(M + l)] + [l + 4p(^ + l)]3/2 



Gn 



24p^ 

-(l + 6p(/i-l)] + [l + 4p(M-l)]-^/" 
24p2 



For weak sixth order nonlinearity p ^ 1, or just very 
close to the bifurcation points, so that ^(/i^ — 1) <C 1 (but 
the effects of dissipation are still small) , the expressions 
( ^ ) , (^) go over into the asymptotic expressions for the 
activation energies (p2|), 5„ cx (p^ — 1)^- 



For strong sixth order nonlinearity, p ^ 1 and p{y? — 
1) <C 1, it is convenient to rescale the dynamical vari- 
ables, 

X - p^'^X, Y = p'/^Y, G = p'l^G (36) 
G = i(p - X)X- + lip + 1)Y^ + y')\ 



and the noise intensity D (|To|), 



8w|,r2 



<j){L0F/2). 



(37) 



The maximal and saddle values of G, Gi,2 and G^, 
respectively, are given by the simple expressions 



Gi,2 = J(P + 1)'/'; 



(38) 



Gu = for - 1< p < 1, Gu = \{p - if'^ for p > 1, 

whereas the expressions for the trajectories in polar co- 
ordinates (i? = p^^^R, (p) are of the form 



^±(G, if) = 2(p-cos2</7)i/2 cos 



OiG,ip)T7r 



, (39) 



9{G,(p) = arccos hc/ip - cos 2(^)^/2 



With account taken of (^), the expressions ( ^ ) for 
the functions M = p^'^^^M, N = p^'^/^N take the form 

M(G) = y"d^[3(p-cos2(^)]i/2sin[6i(G,(^)/3], (40) 

N{G) = - J dip{2p - 3cos2<p) [3{p - cos2(p)]^/^ 
xsin[0(G,^)/3], 



(the above expressions refer to the trajectories which sur- 
round the period two attractors). For the trajectories 
that surround the zero-ampUtude state X ~ Y ^ we 



obtain from (32), (M) 



M(G) 



2tt 



dip [/i — COS '2ip]^^^ ( 



0iG,ip) 



(41) 



27T 



N{G) = - / d^(2/i-3cos2^)[^-cos2(^]i/2 
Jo 

'e{G,ip)+n 



X cos 



6nG. 



The expression for the escape activation energies Sn — 
p^^'^Sn has the same form as Eq. (^8|) for S'„, with 
G, M, N in (|l|) replaced by G, M, N, respectively. The 
escape rate in the variables with tilde has the same form 
as in Eq. (|n|), 



VF„ = Cexp(-5„/^). 



(42) 



It follows from Eqs. (p7|), (p6|), (^) that, in the limit of 
large sixth-order nonlinearity, the fourth-order nonlinear- 
ity parameter 7 drops out of the expressions for the ac- 
tivation energies Sn, the reduced noise intensity D, and 
the escape rates. This is in agreement with Eq. (35), 
which shows explicitly that Sn oc p~^^^ for large p, and 
therefore Sn — p^^'^Sn is independent of 7. 

For large frequency detuning, /i ^ 1, it follows from 
Eqs. (|l|), (|ol), (|4l]) that the activation energies 5i,2,5o 
are of the form 



51,2 ~ Ca^ 



-1/2 



5o~CA^'/^ (/i»l). 



(43) 



It is clear from the asymptotic expressions (^5|), ( }43| ) 
that the activation energy of escape from the period two 
attractors 6*1,2 is a nonmonotonic function of p, i.e. of 
the frequency detuning ujp — 2a;o- The decrease of 61,2 
for large p can be understood as follows. As we men- 
tioned before, the effective reciprocal "temperature" of 
the distribution of the system over the energy G, which 
is given by the M/N = M/N, is determined by the ra- 
tio of the rates of the drift and diffusion of the oscillator 
over G (for a Brownian particle this ratio is indeed equal 
to 1/kT). The drift coefficient is linear in the charac- 
teristic velocity X,Y of the oscillator in rotating frame, 
whereas the diffusion is quadratic in this velocity. There- 
fore when the velocity is large, the ratio M/N becomes 
small. This happens for large p. Here, the characteristic 
time scale for the motion with a given G is set by the 
reciprocal frequency detuning 1/\luf — 2a;o|, and the ra- 
tio M/N oc 1/\lof — 2a;o| oc 1/p. On the other hand, the 
energy interval Gi — G„ for large p is increasing with p 
sublinearly, if the sixth order nonlinearity is dominating 
(and linearly, if the fourth order nonlinearity is dominat- 
ing). Therefore, for large p and p the activation energy 
decreases with the increasing p. 



The position of the maximum of 6*1,2 and the overall 
dependence of Sn on p can be obtained by numerical 
integration of Eqs. (p^), (|4C|), (|4l|). The results are shown 
in Fig. 5. 
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APPENDIX A: WEAK SIXTH ORDER 
NONLINEARITY 

When pC <C 1 the zeroth-order in p( values of the 
functions M(G) and N{G) in the Eqs. ^) and (H) are 
determined by the Eqs.(pO|) and (pT|), respectively. In the 
case of period two attractors to the first order in pC, the 
functions Af , have the form 



-{p ~ cos2(^) 



M(G) « m(°)(G) + A/(^HG),iV(G) « a(")(G) + a(i)(G) 
x[{p-cos2(p)^ -3G], 

n(-^\g)^n[^\g) + n^,^\g), 

n['\g) = -ApC r d^f{G,^)[{p - cos 2^)2 - G] 
J 

fV2 



(Al) 
(A2) 



N 



'\G) = f [2 ((i?f (G,^))^ + {R'^\G,p)f) 

-p ((<)(G, + (i?L"^ (G, ^))5)] //(G, V) 

Here i?^''' are the radii of the orbits with a given G 
evaluated for p = 0. The angles ipi^2 in (|A1[),(|A2[) are 
evaluated in the zeroth order in p(^ 



— i arccos(/i — 2G^/^), (/?2 = tt — 



(A3) 
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FIG. 1. Trajectories and separatrices of the oscillator 
in the absence of noise in slow variables gi , q2 for (a) 
^ = 1.5, n = 0.5 where the stable states of the oscillator 
are period two attractors, and (b) for ^ = 1.5, f2 = 1.5 where 
the steady state q = is also stable. The positions of the 
stable states and the saddle points are denoted by the letters 
tti and s, respectively, (c) The dependence of the dimension- 
less amplitude of the stable (solid line) and unstable (dashed 
line) period two vibrations on O for = 5.0. 
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FIG. 2. Trajectories of the conservative motion (16) for (a) 
fi = 0.5 and (b) fi = 2.0. The values of the Hamiltonian 
function G{X, Y) are shown near the trajectories. The dots 
show the positions of the centers and the hyperbolic points. 
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FIG. 3. The Hamihonian function G{X,Y) (16) (a) for 
/X = 0.5 where the function G has two maxima (which corre- 
spond to period two attractors, with account taken of dissi- 
pation) and a saddle point aX X = Y = Q, and (b) for = 2.0 
where G has two maxima, a minimum at X = F = 0, and two 
saddles which correspond to unstable period two vibrations, 
with account taken of dissipation. 
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FIG. 4. The dependence of the escape activation energy 
Si = 5*2 (lines 1) and So (lines 2) on the scaled frequency 
detuning ^ = Q./C, = 2wf[(wf/2) — lj()\/F in the limit of 
comparatively large fields or small damping, ^ ^ 1. 
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FIG. 5. The dependence of the aetivation energies of the 
phase shp transitions between period two attractors on the 
scaled field amplitude for two different values of the frequency 
detuning f2, as obtained by solving the variational problem 
(11) (solid lines). The dashed lines show the low-damping 
(large C) asymptotes (cf. Fig. 4) 
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FIG. 6. The dependence of the escape activation energies 
Si = §2 (line 1) and So (lines 2) on the scaled frequency 
detuning jj, = Q/( = 2uJF[{tOF /2) — loo]/F in the limit of 
comparatively large fields or small damping, ^ 1, for large 
sixth order nonlinearity, p S> 1. 



